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Abstract. Recent evidence of electronic coherence during energy transfer in 
photosynthctic antenna complexes has reinvigorated the discussion of whether 
coherence and/or entanglement has any practical functionality for these molecular 
systems. Here we investigate quantitative relationships between the quantum yield of a 
light-harvesting complex and the distribution of entanglement among its components. 
Our study focusses on the entanglement yield or average entanglement surviving a 
time scale comparable to the average excitation trapping time. As a prototype system 
we consider the Fenna-Matthews-Olson (FMO) protein of green sulphur bacteria and 
show that there is an inverse relationship between the quantum efficiency and the 
average entanglement between distant donor sites. Our results suggest that long- 
lasting electronic coherence among distant donors might help modulation of the light- 
harvesting function. 



PACS numbers: 71.35.-y, 03.65.Ud, 03.65.Yz 
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1. Introduction 

The first step in photosynthesis is carried out by specialized pigment-protein complexes 
which absorb sunlight and, with near-unit efficiency, transfer the associated electronic 
excitation to reaction centres (RCs) where chemical energy storage begins. Numerous 
investigations have been devoted to elucidate the design principles of natural light- 
harvesting antenna [H [2j [3] and it is currently well known that a key element underlying 
their efficient functioning is the spatial arrangement of their pigment molecules together 
with their electronic interactions. However, the answer to whether electronic coherence 
in these systems could survive, first at physiological temperatures, and second, for long 
enough to influence the energy transfer process, had until recently, remained elusive. 
Remarkably, recent studies are providing evidence that, indeed, some of these fascinating 
molecular aggregates are designed to sustain quantum coherent transfer for longer than 
expected [U [7J [8] , and that they can do so at temperatures as high as in biological 
conditions [HI |7J. A direct implication suggested by these studies is that quantum 
coherence may be affecting light-harvesting processes in vivo [9], but its exact role is 
yet to be understood. 

The observation of coherent excited state dynamics in light-harvesting systems has 
therefore given a new impetus to the long-standing discussion regarding relationships 
between energy transfer efficiency and quantum features such as exciton derealization 
(electronic coherence) and its coherent oscillatory behaviour [T0 | [TT | [12] . In fact, several 
works have discussed highly efficient transport as an optimal interplay of coherent 
dynamics by which exciton phase relationships are preserved, and incoherent processes 
that destroy those relations [HI H31 HH HH HH1 HZl HH HH1 [20]. Within this framework, 
the role proposed for coherent transfer varies from being a strategy to overcome energy 
traps [21] to a mechanism that may speed up the search of the lowest energy level [I]. 

The practical functionality of quantum coherent phenomena may, however, be 
subtler than efficiency maximization. For instance, the structure and arrangement 
of antenna complexes are not solely optimized for fast transfer but also incorporates 
the ability of modulating their function under different environmental conditions 
[22| |23| [2U [251 EE [27J. Then, one can ask how electronic coherence and its possible 
coherent evolution may help the requirement of having built-in mechanisms to modulate 
transfer efficiencies under different environmental parameters. Indeed, it may be that 
light-harvesting systems are able to take full advantage of unique quantum interference 
effects that include the possibility of efficiency regulation according to initial state 
conditions [HI [19] or through accumulative quantum phases in closed transfer pathways 

An outstanding point of the recent experimental results is that the observed 
electronic coherence expands several molecules across the whole antenna complex, 
including distant and weakly interacting pigments [H [6]. This long-range electronic 
coherence has motivated recent studies to investigate electronic energy transfer from 
the perspective of quantum entanglement [281 1291 130] . An intriguing question is 
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whether correlations or interferences between distant pigments play a part in the efficient 
functioning of a light-harvesting complex or whether they are just a consequence of 
exciton derealization with no particular relation to transfer efficiencies. 

In this work we use the framework of entanglement to quantify the strength of 
electronic coherences in a light harvesting system and investigate how such coherences 
are distributed among the molecular sites for different initial conditions and bath 
parameters. Our study focuses on the effects of the coherences that survive a time 
scale that can be of biological relevance such as the average excitation trapping time. 
These long-lasting electronic coherences are quantified by the entanglement yield defined 
as the average value of entanglement at times when trapping events are likely to occur. 
As a prototype system we consider the FMO protein in green sulphur bacteria, where 
coherent dynamics was first observed at 77 K [4] and recently at 277 K [7]. Our 
results suggest that long-lived electronic coherence between well-separated pigments 
may modulate the efficiency profile as a function of the system-bath coupling. This 
modulation is manifested in an inverse relationship between quantum efficiency and the 
entanglement yield among distant pigments acting as excitation donors. We discuss how 
the possible relation between efficiency and donor-donor coherence could be exploited 
in light-harvesting systems. 

2. Energy transfer dynamics in light-harvesting systems 

Under normal operating conditions energy transfer dynamics in light-harvesting 
complexes can be described under the assumption that there is at most a single excitation 
present [251 EH- I n this regime, the general form of the system's density matrix p(t) 
takes the form: 

N N 

p{t) = aoo(*)|0)(0| + J2 a ™n(t)\m)(n\ (1) 

m=l n=l 

where \m) = |0..1 m ....O) denotes the state with the excitation on site m and all other 
chromophores (sites) in the ground state, aoo(^) accounts for any losses that relax the 
system to the ground state |0), and a mn (t) describes the populations (m = n) and 
coherences (m ^ n) in the site basis. 

The coherent evolution of p(t) is dominated by the Hamiltonian describing the 
electronic states of N interacting chromophores (H = 1): 

N N 

H S = e ^ a m°m + Yu V mn(^ a n + a n O ( 2 ) 

m=l n>m 

where e m is the energy of the m — th site, V mn describes the electronic coupling between 
chromophores m and n, and cr+(~) is the usual creation (annihilation) operator of an 
excitation at site m. 

The protein environment, which modulates site energy fluctuations and gives rise 
to decoherence, is usually modelled as a bath of harmonic oscillators linearly coupled 
to the electronic system [32]. All the relevant information of the bath structure and 
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system-bath interaction is then captured by the spectral density J(cj) and its associated 
site reorganization energy, E r = j^°(J(u)/u)duj. However the full details of the 
coupling between the pigments and the surrounding protein environment are not fully 
characterized for these multichromophoric systems. In the FMO, for instance, there 
is an agreement that the reorganization energy values are in the moderate system- 
bath coupling regime with values between E r = 35 cm -1 [33l 134"] and E r = 25 cm -1 
[35] . but less consensus exists on the specific form of the bath spectral density and its 
characteristic correlation time. Reported bath correlation times, span from 35 fs [21] to 
166 fs [33], indicating that the energy transfer can happen under slow bath modulation 
where memory effects may be important. Importantly though, for bath-correlation 
times within this range, long-lived quantum coherence is expected at physiological 
conditions [2~T] . The above discussion supports the idea that an understanding of the 
role of coherence and entanglement under weak system-bath coupling assumptions could 
give an insight into phenomena one can explore in more general conditions for which 
coherence survives. 

Here we consider energy transfer dynamics assuming generators of the Lindblad 
form p{t) = —i[Hs,p(t)] + L(p(t)) which ensures complete positivity. First a 
quantum master equation, derived within the weak-coupling, Born-Markov and secular 
approximations [191 ESI E] is considered: 

Li{ P {t)) =EE7»[4W)4(«) - \{Al(u)A n (u),p}\ (3) 

u) m,n ^ 

This generator accounts for relaxation and dephasing in the eigenstates {|^)} of 
H s , with Lindbland operators given by A m (u) = E^,- £ *=„ <4(*)c m (*') and 
{cu} the frequency spectrum associated to energy differences between single-excitation 
eigenstates. Note we have neglected the Lamb-shift type Hamiltonian. We assume 
that fluctuations at different sites can be correlated with rates given by 7 mn (w) = 
Jo(d mn / \B)lmm(u) [19] . Here Jo is the Bessel function of first kind, d mn the distance 
between sites m and n, A# the bath correlation length, and the single-site rate identical 
for all sites is 7 mm (u;) = = 2ti J {\uj\)\N{—uj)\ with N(u) being the thermal 

occupation number. To compare our results with Ref.[T4"] [37] , we consider an ohmic 
spectral density of the form J(u) = (E r u/uj c )exp(— uj/u c ) with cutoff frequency u c . The 
above model is limited to weak coupling and hence fails to reproduce the localization 
effects and the associated decay in the efficiency for large enough system-bath coupling 
[38]. 

The second model we consider is pure dephasing in the site basis [T6l US], [12] , which 
unlike ([3]), does not take into account the quantum nature of the bath. This model, 
however, describes qualitatively the transition from a purely coherent to incoherent 
regimes of energy transfer as function of the dephasing rate jk- The generator reads 

I*Q>(t)) = Y.lMkP{t)A\ - \{A\A k ,p}\ (4) 

k 

with Lindblad operators in the site basis as A^ = \k){k\. In our calculations, for 
simplicity, the dephasing rate is assumed identical for all sites, i.e. 7^ = "jdeph- 
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The energy transfer dynamics governed by the coherent term H$ and the bath- 
induced irreversible dynamics given by ([3]) or (j3J), which preserve the number of 
excitations in the system, can be interrupted by dissipation to the ground state at any 
of the sites with a rate T m or excitation trapping at the RC with rate K m . Both of these 
processes are described by Lindblad generators of the form (@J with operators = 
and rates 7^ = or 7^ = for dissipation and trapping at the RC respectively. The 
effects of these generators are equivalent to the inclusion of a non-hermitian operator, 
H\ oss = —(i/2) Z)ro = i(r m + Km) <T m cr~, into the system Hamiltonian (T3J 03]. 

3. Quantum Coherence and Entanglement 

Quantum superpositions of the states of a composite system give rise to non-trivial 
quantum features such as coherence and entanglement. These two phenomena are, 
however, not the same. While quantum coherence is generally associated to interference 
effects, entanglement refers to non-classical correlations between distinguishable modes 
or subsystems of a multipartite complex [39] HQ]. However, in scenarios where the 
system's quantum evolution is restricted to a subspace with a single-excitation, quantum 
coherence can be not only a necessary but also a sufficient condition for the existence of 
entanglement [29]. Whether this entanglement can be used as an information resource 
is an ongoing discussion [H], though that is not the concern of this paper. Our purpose 
here is to use entanglement as a framework to quantify electronic coherence and its 
distribution in a light-harvesting complex. 

One of the most striking features of entanglement in the multipartite setting is the 
constrains in the way such quantum correlations can be distributed among subsystems 
[4~2l H3l HI] . This property, known as monogamy, implies that entanglement cannot be 
freely shared among different objects, and therefore, the amount of quantum correlations 
that exist between two subsystems limits strongly the amount of entanglement they can 
share with a third party. The monogamous nature of quantum correlations is expressed 
in strong inequalities bounding the way bipartite entanglement can be distributed, 
as originally demonstrated for pure states of three qubits by Coffman, Kundu, and 
Wootters [12] ■ The generalization of such inequality to the case of pure and mixed 
states of N qubits has also been demonstrated [43] . In particular, it has been shown 
that for a collection of N qubits the distribution of bipartite entanglement, as measured 
by the squared concurrence or tangle [45], satisfies the inequality: 



where r(p n: i...jv) refers to the tangle between the partition of site n and remaining part 
of the system, and r{p mn {i)) is the tangle of the reduced state p mn (t) of two sites m and 
n. 

In the context of the energy transfer dynamics we are interested in, the 
chromophores are considered as effective two-level systems or qubits and the two- 
qubit reduced dynamics p mn (t) is restricted to the zero- and single-excitation subspace. 



N 



E 




(5) 
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Hence, its associated concurrence is simply given by C(p mn (t)) = 2\a mn {t)\ and the 
corresponding tangle by r(p mn (t)) = 4\a mn (t)\ 2 , where a mn (t) are the off-diagonal 
elements of the density matrix in equation [1] representing electronic coherence, showing 
that, in this case, entanglement is a measure of quantum coherence. The inequality 
in equation [5] consequently means that, quantum coherence is also subject to certain 
shareability laws. 

For our purpose of studying quantitative relationships between quantum yield and 
the spatial and temporal distribution of entanglement (or equivalently coherence) among 
the different chromophores in a light-harvesting complex, it therefore appears important 
to consider a global measure that captures the monogamy relation. In order to define 
such a measure one can notice that a summation over all tangles in equation (jSj) leads to 
an inequality for the sum of all pairwise contributions i.e. Ex(t) < \ J2n=i T (Pn) where 

M*)= E r{p mn {t))=A J2 \a mn (t)\ 2 (6) 

m,n>m m,n>m 

Notice that in the case of a pure state, Et^) is closely related to the inverse participation 
ratio, often used as a measure of exciton derealization [IT] . As it has been discussed 
before [281 129], the entanglement of p(t) measured by E T (t) indicates that in the 
case where a single excitation is present in the system, the existence of non-vanishing 
off diagonal elements a mn (t) in equation [TJ is sufficient condition for the presence of 
entanglement. However, we would like to highlight that here the use of this measure 
has been motivated by being one capturing the shareability of bipartite entanglement, 
important for our aim of relating efficiency to distribution of coherence. As a first step 
in understanding such possible relations, we study the entanglement quantified by E(t) 
and its associated yield 0, which will be presented in the next section. 

4. Quantum yield and Entanglement yield 

The performance of a light-harvesting system can be assessed in terms of its efficiency, 
optimality and robustness [3JJ. Here we focus on the quantum efficiency or quantum 
yield 77, which measures the total probability that the excitation is trapped by the 
reaction centre as opposed to being radiatively dissipated at any of the network sites 

ti= / u RC {t)dt (7) 
Jo 

where ujRc(t) = J2 m K m(m\p(t)\m) is the probability density that the excitation is 
transferred from any of the sites to the RC. 

In order to explore quantitative relationships between 77 and the distribution of 
entanglement in space and time, here we concentrate on the amount of entanglement that 
lasts for the time scale of trapping events. In particular, we focus on the expected value 
of entanglement at times when excitation trapping is likely. Formally, this corresponds to 
the average entanglement associated to the waiting time distribution defined by u)Rc{t). 
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Figure 1. (a) Schematic representation of the monomeric subunit in the FMO 
complex. The position of each of the seven individual bacteriochlorophylls (BChls) 
is labeled from 1 to 7. (b-c) Most likely pathways of excitation transfer are shown as 
full arrows: from sites 1 and 2 the excitation is transferred primarily through a direct 
path (b) to the acceptor region, whereas from sites 5 and 6 it is most likely to be 
transferred through an indirect path (c) involving site 7. 



We define this quantity as the entanglement yield: 
1 r°° 

<P=- E{t)u RC {t)dt (8) 

7] JO 

Different from previously used entanglement measures [2SJI2H], the entanglement yield 
quantifies quantum correlations that may be impacting the transport properties on a 
time scale that compares with the average transfer time. In what follows we will in- 
vestigate how this average entanglement gives insight into the modulating effect that 
long-lived electronic coherences may have on the quantum efficiency. 



5. Entanglement yield distribution in FMO 

Contrary to what it could be assumed for antenna systems without symmetry, the 
spatial organization of pigments in a light-harvesting antenna is not random [2H SSj , 
and generally reflects the functionality associated with either individual pigments or a 
set of them. For instance, the FMO complex serves as a link between the chlorosome, 
which harvests the sunlight, and the RC. Theoretical |3Z] and experimental [32] studies 
support the view that each monomeric subunit of the FMO is oriented with sites 1 and 
6 in closer contact with the baseplate chlorosome and with site 3 towards the RC. It is 
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then likely that the excitation enters the FMO by this route, from where it is transferred 
to the target site 3. 

As a starting point to analyze the distribution of the entanglement yield (jHJ) in the 
FMO, we partition the system based on the predicted transfer pathways at 77 K [331H9] . 
The results reported in [33] indicate that there are two preferred non-cascade exciton 
relaxation pathways. One path involves mainly transfer from the strongly interacting 
pair 1 and 2 to the acceptor region defined by the sites 3 and 4. The second pathway 
takes the excitation from the sites 6 and 5 through site 7 until it reaches the acceptor 
region. We therefore define the pairs involving pigments 1-2 and 5-6 as donor dimers, 
whereas the 3-4 pair is denoted the acceptor dimer. The possible transfer pathways 
are classified between direct and indirect transfer pathways, depending on whether they 
exclude or include pigment 7 as illustrated in figure 1(b)- (c). With these partitions in 
mind, we investigate how the total entanglement Exit) distributes along these pathways 
by partitioning Exit) into donors, acceptors, intra-dimers, and mediator contributions: 

E T {t) = E DD + Ead + E dim + E 7trest (9) 

where 

E D D = 7-1,5 + 7"l,6 + 7- 2 ,5 + 7- 2 ,6 (10) 

E A d = T i,3 + T iA ( n ) 

i=l,2,5,6 

Edim = 7-1,2 + 7-3,4 + 75,6 (12) 

Eldest = ^2 T 7,i (13) 
i+1 

Notice that the inter-dimer entanglements donor-donor Edd an d acceptor-donor Ead, 
account primarily for long-range coherences while the mediator entanglement Ej trest 
and the intra-dimer entanglement E^m quantify mainly short-range coherences. Their 
corresponding entanglement yields are defined as 4>dd, (ft ad, 4>dim and (p-j^est and are 
calculated through (jHJ). Note that for ease of notation r m , n = r(p mn ) in (I10lll3p . 



In ([3]) we assume the same parameters as in [2J with a thermal bath at room 
temperature T = 293i^ and a cutoff frequency u c = 150 cm -1 . The excitation is 
assumed to be transferred to the RC from site 3 with a rate K3 = 1 ps" 1 and radiative 
relaxation occurs with a rate of T m = V — 1 ns _1 , identical for all pigments [31]. The 
site energies and electronic couplings in (J2j) have been taken from table 1 in reference 
[33]. 



5.1. Entanglement survival 

In order to illustrate the entanglement dynamics, we plot in figure [2]the time evolution of 
the different entanglement contributions i.e. E D d, Ead, E dim and E 7)J . est , assuming the 
excitation initially localized on site 6 and decoherence described by model (J3j). Here we 
consider E r = 35 cm" 1 and uncorrelated energy fluctuations (As = 0). For comparison, 
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Figure 2. Weak- coupling and secular approximation model. Time evolution of the 
contributions to the total entanglement (-) and RC trapping probability density (- 
-) for excitation initially localized on BChl 6, at T = 293 K, E r = 35 cm" 1 and 
A b = 0. (a) Edd, distant-donors contribution (b) E^im, intra-dimer contribution (c) 
Ead, acceptor-donor contribution and (d) E 7trest , site 7 contribution. 



the trapping probability density U)Rc(t) (dashed line) is also shown. Several interesting 
features can be drawn from figure [2j Firstly, the total entanglement, which is the 
sum of these four contributions, is predicted to have a maximum value first at around 
0.04 ps and then at around 0.8 ps, close to where the trapping probability reaches a 
maximum too. These time scales of non-vanishing entanglement compare well with the 
experimental observations of coherence surviving at least 0.3 ps at 277 K [TJ. Secondly, 
the entanglement values that we are interested in correspond to the overlap with the 
waiting time distribution ojRc{t), that is, the entanglement at times where there is non- 
vanishing trapping probabilities, i.e. between t = 0.01 ps and t = 2 ps (see figure |2]). 
Notice that although the entanglement values are small, there is a considerable overlap 
in time where both ujRc{t) and the entanglement contributions have values different 
from zero. Interestingly, except for Edd, the entanglements Ead, E 7)rest and E^ m peak 
closely where ujRcit) does. With the exception of the second peak in E dim , these results 
are similar to the case where the excitation is initially localized on site 1 (results not 
shown). Given that the quantum yield is near unity with the parameters considered 
(figure Efa)), we can conjecture that larger efficiencies are associated with both lower 
yet not zero values of Edd and larger values of Ead and E 7 ^ rest at typical trapping times. 
In order to explore this hypothesis, we now turn to discuss the entanglement yield. 
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5.2. Efficiency and donor-donor entanglement yield 

Quantitative relations between efficiency and the different contributions to the total 
entanglement yield were investigated by comparing their values under different bath 
parameters and various initial conditions. The initial states considered are excitation 
localized on either site 1 or 6. Under decoherence given by (0) with uncorrelated 
fluctuations (As = 0), these two initial states exhibit, in general, different efficiencies as 
a function of the system-bath coupling quantified by E r (see figure 0^a)). As previously 
reported [H], this model predicts a monotonic increase of r] with increasing E r , which is 
a behaviour valid only for weak system-bath couplings [38]. Therefore, the maximum E r 
we consider is E r = 50 cm -1 . As interaction with the environment is stronger, coherence 
is destroyed and, consequently, the total entanglement yield decreases monotonically 
as shown in figure 0^b). Interestingly, the interplay between coherent and incoherent 
processes prompts a crossing of the efficiency profiles for the two different states at 
around E r « 0.021 cm -1 , indicating that one state is more efficient than the other 
one depending on the coupling to the environment. A similar crossing is observed 
for the total entanglement yield but at a slightly larger value i.e. E r « 0.035 cm -1 . 
However, such crossing is not necessarily exhibited by all the contributions to the total 
entanglement yield as it can be seen in figure 0(c)- (f). In particular, notice that 4>dd 
shown in figure 0(c) is the only contribution featuring a profile crossing at around 
E r « 0.017 cm -1 . The most important feature of figure0is the fact that for a large range 
of E r parameters an inverse relationship between r] and 4>dd is suggested. For instance, 
if we compare the efficiency and entanglement yield values for the two initial states at 
E r = 10 cm -1 , we can see that they satisfy ri l {E r = 10cm -1 ) < r] 6 (E r = 10cm -1 ) 
while 0^ D (-E r = 10cm -1 ) > <fP DD (E r = 10cm -1 ). This supports the idea that higher 
efficiencies are associated to lower values in donor-donor entanglement. 

Importantly, the trends calculated for <fidim, 4>ad and 4>7, r est as function of bath 
reorganization energy, E r (see figure 0(d)- (f)) illustrate the differences in the distribution 
of entanglement among sites depending on the initial state. For example, the average 
entanglement of mediator site 7 and the rest, <pi^esti is always larger for the case when 
the excitation is initially localized on site 6 than for the initial state at site 1 since, as 
discussed above, energy migration from site 6 is largely mediated by site 7 (see figure 
1(c)). Similarly, the donor-acceptor entanglement yield, <Pad is in general larger for 
initial state at site 6 than 1 as interactions between sites 5 and 4 are larger than between 
sites 1 or 2 and sites 3 or 4 (cf. hamiltonian given in [33]). In contrast, <pdim is always 
greater for the initial state at site 1 consistent with the fact that energy transfer from 
this site has a large contribution of the coherence associated to the strong interacting 
1-2 pair. We would like to point out that so far we have not found clear relationships 
between efficiency and the contributions (f)di m , 4>ad and 4>i,rest- 

The majority of the above analyses and, in particular, the identified inverse 
relationship between §dd and rj hold when decoherence is given by local random 
fluctuations described by the pure dephasing model (jlj) as illustrated in figure HI The 
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E [cm '] E [cm '] 

r r 



Figure 3. Weak- coupling and secular approximation model, (a) 77, quantum yield (b) 
<f>T, total entanglement yield (c) (j>DD, donor-donor entanglement yield (d) <pdim, intra- 
dimcr entanglement yield (e) 4>ad, acceptor-donor entanglement yield and (f) <f>j t resti 
site 7 entanglement yield, as a function of reorganization energy, for the initial states 
localized on site 1 (-) and site 6 (- -), for T = 293 K and As = 0. 
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quantum yield and the total entanglement yield as functions of the dephasing rate jdeph 
are shown in figures H(a) and IU(b) respectively, while the different contributions to 
the entanglement yield are depicted in figures H^c)-(f). As reported in [161 [IS], the 
quantum yield has a non-monotonic behaviour as the dephasing rate increases (see 
figure H(a)), being qualitatively consistent with a transition from coherent to incoherent 
transfer. The efficiency profiles for the two initial states cross (figure S(a)), similarly to 
the weak-coupling results (figure [3](a)). Nevertheless, there is no intersection in their 
total entanglement yields (compare figures EJ^b) and|4(b)). Entanglement yield crossing 
is however observed in the donor-donor contribution (compare figures [3]^ c) and|4](c)), 
pointing out again the inverse relationship between quantum efficiency and the average 
of donor-donor entanglement at trapping times. 

The main difference between the two decoherence models here considered is the 
crossing of the profiles of the total entanglement yield for the two initial states, which 
for the secular case is observed where coherence has not vanished (see figure [3]^a)), while 
in the pure dephasing model there is no such crossing before coherence is practically 
zero (see figure H](c)). We believe these differences arise from the bath features that 
are considered by each model. It has been discussed that the phonon-assisted transport 
properties in FMO complex are greatly dominated by the temperature-independent 
spontaneous emission of energy into the phonon bath [TJ]. These fluctuations are 
accounted for within the secular approximation while they are completely neglected by 
the classical bath approximation in pure-dephasing models. Despite of the differences, 
however, the suggested inverse relationship between the efficiency and the average of 
donor-donor entanglement 4>dd holds for the two decoherence models. 

5. 3. Entanglement yield under spatially correlated fluctuations 

Long-range spatial correlations have been suggested as an important mechanism helping 
quantum coherence to survive at high-temperatures in light-harvesting systems (H El [7J . 
This hypothesis is supported by recent studies showing that for a single donor- accept or 
pair environmental fluctuations that are correlated in space lead to a high critical 
temperature for which a transition from coherent to incoherent transfer occurs [50] . 
In this section we explore the robustness of the above results when bath-induced energy 
fluctuations of different pigments are correlated [19]. In particular we are interested in 
determining whether the suggested inverse relationship between 77 and 4>dd still holds. 
We therefore fix the reorganization energy to E r = 35 cm -1 and compare the profiles 
of quantum and entanglement yields as functions of the degree of bath correlation by 
varying A# in ([3]). In the presence of spatial correlations the efficiency for the initial 
states considered decreases [37] and distinguishability of the initial state becomes evident 
[T9] (see figure [5]( a)). As correlations are included the coherent nature of the transfer 
dominates and the entanglement yield increases. Importantly, the inverse relationship 
between efficiency and the average entanglement among donors is again observed (see 
figure [5](b)). 
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Figure 4. Pure dephasing model, (a) 77, quantum yield (b) 4>t, total entanglement 
yield (c) 4>DDi donor-donor entanglement yield (d) 4>dim, intra-dimer entanglement 
yield (c) 4>ad, acceptor-donor entanglement yield and (f) (pT^est, site 7 entanglement 
yield, as a function of dephasing rate for the initial states localized on site 1 (-) and 
site 6 (- -). 
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Figure 5. Weak- coupling and secular approximation model, including spatial 
correlations, (a) Quantum yield and (b) 4>ddi donor-donor entanglement yield as 
a function of bath correlation length, for the initial states localized on site 1 (-) and 
site 6 (- -), for T = 293 K and E r = 35 cm" 1 . 



6. Concluding Remarks 

A key challenge in the theory of photosynthetic energy transfer is formulating possible 
quantitative relations between electronic coherence and the light-harvesting function. 
In this work we have approached this challenge by investigating relationships between 
the transfer efficiency and the entanglement present in the system at times when 
trapping events are likely to occur. We have shown that in the FMO complex and 
under weak system-bath coupling assumptions, this long-lived coherence is spatially 
distributed in such a way that the average entanglement among donor dimers exhibits 
an inverse relationship with the quantum yield. This feature is observed under different 
decoherence models and a wide range of bath parameters including room temperature, 
suggesting that donor-donor correlations may be exploited for efficiency modulation. 
Further investigation is however needed in order to asses the validity of this quantitative 
relation in more general theories of energy transfer that are not limited to the weak- 
coupling limit here explored [381 EB ED E2] . 

The ability to regulate transport properties under different environmental 
conditions is crucial for robust photosynthetic systems. Under high light-conditions, 
for instance, photosynthetic systems switch on photo-protection mechanisms that allow 
them to control (reduce) how fast excitation energy is converted into chemical energy 
[22l [23] . The mechanism here described by which a convenient distribution of electronic 
coherence among pigments could allow control of the efficiency profile, may prove useful 
in adjusting transport properties under variations of light-intensity. It will therefore 
be interesting to investigate correlations between efficiency and the spatio-temporal 
distribution of coherence when pigment-protein complexes are exposed to different light 
intensities. 
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